Critical exponents of the degenerate Hubbard model 



Holger Frahm* 

Institut fur Theoretische Physik, Universitat Hannover 
D-3000 Hannover 1, F. R. Germany 

Andreas Schadschneider^ + 

Institut fiir Theoretische Physik, Universitat zu Koln 
D-5000 Kdln 41, F. R. Germany 

July 15, 1992 

Submitted to J. Phys. A 
Abstract 

We study the critical behaviour of the SU(iV) generalization of the one-dimensional 
Hubbard model with arbitrary degeneracy N. Using the integrability of this model by 
Bethe Ansatz we are able to compute the spectrum of the low-lying excitations in a large 
but finite box for arbitrary values of the electron density and of the Coulomb interaction. 
This information is used to determine the asymptotic behaviour of correlation functions 
at zero temperature in the presence of external fields lifting the degeneracy. The critical 
exponents depend on the system parameters through a, N x N dressed charge matrix 
implying the relevance of the interaction of charge- and spin-density waves. 
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1 Introduction 



The physics of highly correlated electron systems has long been the subject of extensive studies 
in condensed matter physics. Recently, the non-Fermi liquid character of low- dimensional sys- 
tems has attracted renewed interest in one- dimensional realizations of these systems where large 
quantum fluctuations lead to Luttinger liquid behaviour [0J: the correlation functions decay as 
power-laws at zero temperature, the exponents depending on the system parameters such as 
electron density, magnetization (or applied magnetic field) and strength of the interaction. 

In this context exactly soluble models can provide a variety of new insights, in particular 
when used together with the general results on quantum critical behaviour in one spatial di- 
mension as provided by the theory of conformal invariance @-@]. Here the universality class of 
the quantum system is completely determined by a single dimensionless number — the so-called 
central charge c of the underlying Virasoro algebra. This number c as well as the dimensions of 
the operators present in the theory can be extracted from analytical results for the spectrum of 
low-lying states in finite geometries. In this language, Luttinger liquids correspond to a central 
charge c = 1, the dependence of the anomalous dimensions on the system parameters is through 
a single dimensionless number — the coupling constant of the corresponding Gaussian model. 

The situation described above is the generic behaviour of ID quantum systems with a 
single critical degree of freedom, as realized in spin-| chains or systems of spinless fermions. 
The situation becomes more complicated for systems where two or more massless excitations 
are possible: the lack of Lorentz invariance (the corresponding Fermi velocities differ in general) 
prevents the direct application of the predictions of conformal field theory and the interacting 
nature of the system complicates the factorization of the problem into independent ones for 
each critical quasi-particle mode. On the other hand the exact results on the finite-size scaling 
of the low-lying energies available for Bethe Ansatz soluble models suggests a resolution of this 
problem: the spectrum is that of a mult i- component Gaussian model. In analogy to the analysis 
of conformal invariant theories the universality class of a system with N massless collective 
excitation modes is determined by N dimensionless numbers c r — reducing to the central charge 
in the scalar case. The anomalous dimensions are functions of the system parameters not 
through a single coupling constant but through & N x N matrix of dimensionless numbers — 
the so-called dressed charge matrix. This behaviour has been found in a large number of ID 
quantum systems, including certain integrable spin-chains with S > | |J || , and the Hubbard 
J7]. H, [| and t— J-models of correlated electrons, the latter both at the integrable supersymmetric 
point J/t = 2 [ID| and away from integrability |]TT . 



Most of the integrable models in this list are solved by a hierarchy of Bethe Ansatze fll2 
The first one introduces a set of wavenumbers describing the phase of the wave function and 
determining the spectrum, the others are necessary for the wavefunctions to show the symmetry 
corresponding to a particular representation of the Permutation group. This allows for the 
solution of certain systems with various choices of internal degrees of freedom, e.g. an SU(2)- 
spin in the Hubbard model. 

In the present paper we study the critical properties of a generalization of the Hubbard 
model, describing electrons carrying an SU(iV)-spin index on the lattice. Unlike the case of the 
regular (N = 2) Hubbard model this model allows to study a Mott transition at a finite value 



of the Coulomb interaction |TJ, |15f . Our paper is organized as follows: In the following 
section we shall introduce the model and present a brief review of its Bethe Ansatz solution and 
a qualitative discussion of its excitation spectrum. In Section 3 the results for the finite-size 
corrections to the energies of low-lying states as well as their relation to the critical exponents 



are given in terms of the NxN dressed charge matrix. In Section 4 the integral equations for this 
matrix are solved in the zero-field case and the operator dimensions are computed as a function 
of the electron density and the strength of the interaction. They are shown to reflect the full 
SU(iV)-spin symmetry present in this case. These results are applied to the computation of the 
critical exponents for some correlation functions of interest. Finally, in Section 5 we consider 
states where the SU(iV)-symmetry of the ground state is broken by magnetic fields coupling to 
the various flavours of the internal degree of freedom. In the limiting case of strong coupling 
we discuss the dependence of the dressed charge matrix and the critical exponents on these 
external fields. 



2 The Bethe-Ansatz solution of the degenerate Hubbard 
model 

The Hamiltonian of the degenerate Hubbard model on a chain of length L is given by the 
following expression: 



L N L 

n N = - Y,H V { c )+i,sCj S + c] s c i+ i, s ) v + 4m X n n 

j=l s=l j=l s,s'(sj^s') 

L N L N 

+ t*J212 n is + Y,Y, h s n 3s (!) 

j=l s=l j=l s=l 

The Fermi operator cj s (cj S ) creates (annihilates) an electron at site j with spin index s G 
{1,...,N} and rij s = c^ s Cj S is the corresponding number operator. The real parameters h s 
may be considered as generalized magnetic fields. Note that these fields do not destroy the 
integrability of the model since the numbers N s = Ylf=i n js of particles with spin index s 
are conserved. V projects onto the subspace of states having at most two electrons at each 
site. This projection is necessary in order to render the Hamiltonian ([!]) Bethe-Ansatz solvable 
[0], [T7|, (TP| and has no effect for N = 2. The Bethe-Ansatz wave function which solves the 



Schrodinger equation for ([[]) for a total number iV c of electrons is characterized by the momenta 
kj (j = 1, . . . , N c ) and N — 1 sets of rapidities (s = 1, . . . , N — 1; a = 1, . . . , M s ). Imposing 
periodic boundary conditions on the wave function leads to the Bethe-Ansatz equations 
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Here we have set M = iV c , Mjy = and = sin kj. 

The quantum numbers Ij and are integer or half-integer depending on the parity of the 
numbers N c , 

T M l T(s) M s ~ M s-1 ~ M s+1 + 1 

Ij = — mod 1, J [ a s> = mod 1. (3) 



Energy and momentum of the model in a state corresponding to a solution of (0) are completely 
determined by the momenta kf. 

N c N 

E = -2 cos kj + pN c - Y h s N s, 

3=1 s=l 
N c 9 / N c N-l M s \ 

p = E*i = r E^ + EE4 S) , (4) 
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where N s denotes the total number of electrons with orbital index s. In the thermodynamic 
limit (L — > oo, with N c /L, M s /L kept constant) the eqs. (0) corresponding to the ground state 
of (HD can be transformed into a set of coupled integral equations for the densities p c (k) and 
p s (X) of the parameters kj and X^\ respectively: 

Pc(k) = + / 1 d\K 1 (smk-X)p 1 (X), 

Z7T 271 J—Ai 

1 rh) 1 Mi 

PiW = 7T dkK x (X - smk)p c (k) - — / dp K 2 (X - p)p 1 (p) 

1 /" A 2 

+— / dpK x {X- p)p 2 (p) , (5) 

1 A-l 1 /"As 

P»W = 7^-/ dpKi(X- p)p s ^i(p) - — dpK 2 (X- p)p s (p) 

2,71 ,/-A s _i Z7T J— A a 

+ [ Aa+1 dpK^X- p) Ps+1 (p), (s = 2,...,iV-l) 

•>-A s+ i 

with A at = 0. The kernels 1^2(2;) of these eqs. are given by 

The values of the parameters fco and Ai, . . . , Ajy-i are determined through the normalizations 

rko 

n c = dkp c (k), 



-fed 

n s = / As_1 dX p a _ x (A) - dX p s (X) , (s = l,...,iV-l) (7) 

where n c = N c /L is the total density of electrons and n s = N s /L = (M a _i — M s )/L is the 
density of electrons with index s. Furthermore we have set p = Pc and A = fco-Q 
The ground state energy per lattice site is 

rk ^-1 m s 

eoo = / dk(p + hx -2 cos k)p c (k) + 22 (h s+ i - h s ) dX p s (X) (8) 

feo s= i •'—A s 

which may alternatively be expressed in terms of the dressed energy 



, , dfce c (A;). (9) 



1 In the following we will adopt the convention that an index '0' stands for 'c'. 



Here e c (k) is the solution of the system of coupled integral equations 



e c (k) 
s 1 (X) 

e s (X) 



4 0) (^) + t^ d\Ki(wnk-X)e 1 (X), 
2n J-Ai 



Ex (A) H / dk cos k Ki(X — sin k)e c (k) / rf/ii^ 2 (A — jj)si(jj) 



1 f A 2 



2lT J-A 2 



2/T J-fc„ 

djjiKi(X - /i)£ 2 (aO, 



27T J-Ai 



(10) 



e f\^) + 7rl dfiKi(X- //)e s _i(//) - -!- / diiK 2 {X - h)e s {h) 

27T J-As-! 271 J-A s 



1 /"A s +i 

+ / dfiK^X - ii)£ s+ i(h) 
2n J-A s+1 



N-l). 



The bare energies are from 



e^'(k) = n + /ii - 2 cos fc, 4 Uj (A) = /i s+1 - h s . 
The dressed energies ( |10"D obey the conditions 

e c {k ) = 0, e s {A s ) = 0. 



(11) 



(12) 



The ground state at half-filling (n c = 1) and with vanishing fields h s shows for N > 2 
an interesting behavior which has not been found in the case of the standard Hubbard model 
(N = 2). This has been noticed independently by Schlottmann \TA\ and one of the authors 
||15|| . For u > u c one finds ko = tt and for u < u c one has ko < 7T. Here the critical value u c is 
determined through the implicit equation 
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where Gn(x; u) in terms of the Digamma-function if)(x) is given by 
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For iV = 2 we have u c = (see Fig. 1) as already shown by Lieb and Wu [21 



The excitation spectrum in zero-fields has also been studied in []T5| using an extension of 
the method developed in |23| and in [^4|. One finds N — 1 gapless spin excitations with 
soft modes with wave numbers 2sPf (s = 1, . . . , N — 1, Pp = j^n c ) and so-called particle-hole 
excitations.^ These also are gapless and do exist only for k < tt, i. e. for n r < 1 or n r = 1 



and u < u c . In |l5j also excitations corresponding to doubly occupied sites have been studied. 
These are described by complex momenta k^ satisfying sin /c ± = A ± iu and have a finite gap 
(at least for n c — 1 and u > u c ). 

The special structure of the ground state for n c = 1 leads to interesting properties of the 
model. For u > u c the only possible charge carrying excitations are those involving complex 
momenta. As these excitations have a gap the system is in an insulating state. For u < u c 
particle-hole excitations become possible. These excitations may carry a current and so the 
system is in a metallic phase. This shows the existence of Mott-transition at the critical value 

2 For a discussion of the zero-field excitation spectrum for N — 2, see e. g. p2| and references therein. 



u c of the Coulomb repulsion u. The transition is also reflected in the behavior of other physical 
quantities, e. g. the charge susceptibility Xc and the Fermi velocity v c 
excitations with double occupations does not vanish in the limit u \ u c 
to the charge carrying shows a discontinuity at u = u c . We may thus say that the transition is 
of 'first order'. 

In the limit u 
chain 



14] , [TJj . As the gap to the 
the value of the gap 



oo at n c = 1 the model ([[]) becomes equivalent to the SU(./V) Heisenberg 
24]. This equivalence generalizes the well-known relation between the regular (N = 2) 



Hubbard model and the Heisenberg antiferromagnet. 

Recently, Schlottmann |20| showed that in the continuum limit the particles interact via a 
potential of the form 1/ sinh r where r is the distance between the particles involved in some 
properly chosen units. This reflects the nonlocal character of the interaction as introduced by 
the projectors V in the kinetic terms of ([!]) for N > 2. 



3 Finite-size corrections and conformal properties 

As shown in the preceding section the degenerate Hubbard model supports gapless excitations in 
general. Thus we may apply the concepts of conformal field theory to determine the asymptotic 
behavior of the correlation functions, e. g. the critical exponents. 

First we calculate exactly the finite-size corrections to the ground state energy and the 
energies of the excited states. This can be done by a straightforward extension of the calculation 
for the case iV = 2 @. The results can be expressed in terms of the iV x iV-dressed charge 

/ £cc{ko) &i(Ai) ••■ £ c ,WAjv-i) \ 
£ic(fco) fn(Ai) ••• Ci,n-i(An-i) 

: : : • ^ 15 > 

\ 6v-i,c(&o) 6v-i,i(Ai) ••■ 6v-i,Ar-i(Ajv-i) / 
The elements of Z can be obtained from the dressed-charge functions £ rs (A) which obey a 
system of coupled integral equations similar to (|10l). For r = c, 1, . . . , N — 1 we have 
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The finite-size scaling behavior of the ground-state energy is found to be 

E - Le^ = -— v s 



,N-1). 
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with the Fermi velocities of charge and spin excitations 

1 1 



v = v c 
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Energies and momenta of the excitations scale as 
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V = diag(u c ,t>i, . . . ,v N -i). 
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+ 2££P r P F , s+1 , (19) 
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Here N^, Nf are positive integers and AM and D are vectors characterizing the excited state 
under consideration. Pp s are the Fermi momenta for electrons with spin index s. For the 
ground state in the thermodynamic limit we have AM S = 0, D s = (s = c, 1, . . . , N — 1). 
For an excited state AM has integer components denoting the change of the total number of 
electrons and the number of electrons with index s with respect to the ground state. D s are 
integer or half-odd integer depending on the parities of the AM S . Due to ([|) we have 

AN C + AMi 
D c = modi, 

AM S _! + AM s+l 

D s = modi (s = l,...,JV-l) (20) 

with AM = AN C and AM N = 0. 

In general, all velocities v s are different. In this case the results (0) and ([ID may be 
interpreted in terms of a semidirect product of N independent Virasoro algebras.0 All these 
Virasoro algebras have central charge c s — 1. For vanishing fields h s all magnon velocities 
Vi, . . . ,vn-i are equal [H], [T5| and we have a semidirect product of a c = 1 Gaussian theory 
- reflecting the U(l) symmetry of the charge sector - and a c = N — 1 Wess-Zumino-Witten 
theory - reflecting the SU(iV)-symmetry of the spin sector. 

Comparing (|T9|) with the predictions of the conformal field theory |3|. [|] 

9 N-l 

£(AM,D)-£ = -J2^(At + A;), 

L s=0 

cs N-l N-l s 

P(AM,D)-P = ^(A+-A;)+2$:^P r P F , s+1 , (21) 

^ s=0 s=0 r=0 

one obtains expressions for the conformal dimensions Af of the primary fields in terms of the 
dressed charge matrix. Requiring that all dimensions are positive we find 

1 \ 2 

(Z^D) ± - (Z^AM) ) +2^ (22) 



2Af- 

which, in general, depend on the system parameters. In the following section we will show that 
for vanishing fields the Af (s = 1, . . . , N — 1) are functions of the components of AM and D 
only, whereas A^ depends on the strength u of the Coulomb repulsion and the density n c of 
electrons. 

We now make further use of the results of conformal field theory to write down the correlation 
functions for primary fields as 

(0 A ±(a;,t)0 A ±(O,O)) = [[ . + — . ■ (23) 



(x — iv s t) 2As (x + iv s t) 



3 See B and references therein for a more detailed discussion of this point. 



The correlation functions of the physical fields consist of a sum of terms (p3|). In the following 
we will study correlators of the form (Oj(t)Ol(0)) where O is given in terms of c and c^. To find 
the asymptotic behavior of the correlator one has to expand O in terms of the conformal fields. 
This is not possible in general, but the explicit form of O allows for an identification of the 
quantum numbers M c , Mi, . . . , M^-i of the intermediate states. Therefore the leading term in 
the asymptotic expansion of (Oj(t)Ol(0)) can be obtained from (|22|) through minimizing with 
respect to the D s satisfying fl20|) . 



4 Critical exponents for vanishing fields 

In the absence of magnetic fields it is easily seen that A r = oo for all r. This allows to eliminate 
the A-dependent quantities from the Bethe Ansatz integral equations by Fourier transformation. 
From Eq. (|T6|) we obtain for the fc-dependent entries of the dressed charge matrix 



N — r 1 r z o 
£rc(z) = AT + — / dy£ rc (y)G N (z - y; 1). 
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(24) 



where we have introduced a new variable z = smk/u. 

The solution for z c = £ cc (zo) is obtained by iteration for small values of z where we obtain 



Z r 



7 + i>a) 



Zq 
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For Zq 3> N a perturbative scheme |25[ based on the Wiener-Hopf method can be applied, 
giving 

z r ~ VN (l - ^r^-) for z > N. (26) 



2ttzq 

For intermediate values of z$ the integral equation (pl[) is easily solved numerically. The de- 
pendence of z c on the density and the strength of the Coulomb interaction is shown in Fig. |] 
for some values of N. In addition to the interpretation of the z c as measure of the reordering 
of the Fermi-sea due to the interaction when an electron is added there exists a direct relation 
to physical observables: It can be expressed as 

z 2 c = nv c n 2 c K (27) 

in terms of the compressibility k = —(l/L)dL/dp of the electron gas (p being the pressure). 

As in the regular Hubbard model (N = 2) one can employ the Wiener-Hopf method to 
compute the remaining elements of the dressed charge matrix ( |i~5"t ) yielding: 
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Note that the symmetric (N — l) x (A^ — l)-block of this matrix is completely determined by 
the SU(A)-spin-symmetry of the system: this symmetry is manifest in the kernel of the integral 



equations and allows for the reduction of the corresponding matrix Wiener-Hopf problem to 
N — 1 scalar ones [26, p^] . The latter are soluble by quadratures and one obtains closed 
expressions for the matrix elements of Zjj : 



[ N )» ~ N Sm 2N 4 + zl-2 cos ^z r z s - sin 2 g ' ^« " cos ' 
The square of i^ 1 is the Cartan matrix for the Lie algebra SU(iV): 

(-2w 2 ) rs = (Civ)rs = 2d~ r ,s — ^r,s+l — <Wl,s- 

Using the properties of the dressed charge matrix we find for the critical exponents 

/N ~ l A /. iV-S 2 
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For vanishing magnetic field the magnon velocities v r (r = 1, . . . , N — 1) are identical, hence 
only the quantities 

N-l i iV-1 

2A± = 2 J] A± = -AN^AN, + D^D CT ± AN^ • D CT + 2 ]T 7V±, 

r=l 4 r=l 

TV — r 

where (AN ff ) r = AM r - -^AiV c , (D CT ) r = £> r , r = 1, . . . , N - 1 (31) 

appear in the exponents describing the asymptotic behaviour of the correlation functions. Note 
that they are independent of the quantity z c incorporating the dependence of the anomalous 
dimensions on electron density and strength of the interaction. They are completely determined 
by the SU(iV)-symmetry of the zero field ground state and, in fact, of the same form as the 



exponents characterizing the SU(iV)-symmetric critical vertex models and spin-chains ]26, ^7 
(the difference being the possibility of fractional values for the elements of AN CT ). 

Now we are able to study the asymptotic behaviour of correlation functions of interest: 
For the field-field correlation function we have AN C = 1 and AM = (1, . . . , 1*., 0, . . .). The 
corresponding values of D can be read off from Eq. (^0|), the contribution at wavenumber 
k = Pp = (-7T / N)n c arises from the choice D r = \{5 rs — 5 r , a +i) for any s giving 



Hence, the singularity of the momentum distribution function at the Fermi-point 

(c k>s cl,s) ~ \k ~ P F \ a (33) 

is characterized by the exponent 

6 1 1 

a = W? + ~e~N (34) 



where 6 = 2z 2 varies between 2 and 2N as u decreases from oo to 0. Hence, 



= a (u -> 0) < a < a(u -> oo) = - - + — ^ (35) 

showing the Luttinger liquid character of the degenerate Hubbard model at nonzero interaction. 
Numerical data for the exponent a are presented in Fig. |3|. 

For the density-density correlation function (n s (x, t)n s (0)) — (n) 2 we find contributions at 
wavenumbers k = 2mPp (m = 1, . . . , N) with exponents 

„ A+ ^ f m \ 2 „ „ A + m(N — m) . . 

2A " = 2 [m) e > 2A " = (36) 

arising from the choice AN C = AM r = and D r = 5 r ^N-m- In addition, there are k = 0-terms 
decaying as x~ 2 asymptotically. They are generated by the marginally relevant secondary 
operators in the conformal family of the unit operator (i.e. AiV c = AM r = and D c = D r = 
but or one of the in Eq. fl3"0| ) equal to 1). 

From the discussion in the preceeding section it is clear that the above statements are valid 
for any value of the coupling constant and for any filling with n c < 1. For n c = 1 two cases 
have to be discussed separately: 

(i) For n c = 1 and u < u c the system is in a metallic phase (described by k$ = k (u) < n 
in the Bethe Ansatz equations). The critical exponents of the system are given by fl30D with z c 
being a function of the Coulomb interaction through (0). The number z c decreases from \H\l 
to 1 as the strength of the Coulomb interaction is varied from to u c . It is possible to expand 
the integral equations for the density p c (k) in the neighbourhood of the Mott transition u < u c 
to determine ko(u) from <^). For iV — > oo the resulting expressions simplify, giving u c = 



131 and 



z c «l + %^/l-- + ••• (37) 




The same square-root singularity in the Coulomb coupling u — but with a different numerical 
prefactor — is found near u = u c for finite N. Through (|30D it also appears in the critical 
exponents near the Mott transition, e.g. the exponent a of the momentum distribution (|34| ) 
varies like 

a«--4-Jl-- +... (38) 



as u approaches u c from below (for iV — > oo). 

(ii) For u > u c with one particle per site the system is in an insulating state, hence charge 
carrying excitations develop a gap. Excitations in the spin degrees of freedom, however, continue 
to be massless at zero temperature. The critical properties of this state can be described along 
the lines of the discussion above: For n c = 1 and u > w c we have ko = it. Hence, the fc-dependent 
quantities disappear from Eq. (|lfj) leaving a system of iV — 1 coupled integral equations for the 
spin components of the dressed charge matrix. The expression for the conformal dimensions 
are of the form (p2~D and depend on the applied magnetic fields. The correlation functions for 
states corresponding to critical excitations are given by (p3|) . These are exactly those found 
in the SU(iV)-generalization of the Heisenberg spin chain (see e.g. 0). Note, that ( p0| ) implies 
that the momentum of the intermediate state is shifted by 7r for states with odd AMi. 



5 Magnetic field effects in the strong coupling regime 



Nonzero magnetic fields h s in flTT|) lead to finite values of the parameters A s through (|12|) . This 
effect in turn leads to a general dependence of the elements of the dressed charge matrix on 
the system parameters u, n c and all of the fields. Only at and beyond certain critical values 
of the fields where one or more bands are completely depleted the integral equations simplify 
to some extent so that analytical results may become available (for a discussion of the N = 2 
case see 0). A more detailed study of the magnetic field dependence of the dressed charge 
matrix becomes possible in the strong coupling limit ||. As in the SU(2) Hubbard model the 
Bethe Ansatz integral equations (|5|), (|H]) and (|T6| ) describing the system simplify in the limit 
u — > oo (see also In this limit the fc-dependent quantities can be eliminated which allows 

to study the effect of magnetic fields on the critical exponents in more detail. Upon rescaling 
of the variables X/u — > A = A one obtains to leading order 

e r (A) = (A) - f Ar dpK 2 (X - p)e r {p) 

271 J-A r 

+— / dfiKi(\ - fi)e r -x(fJ>) + 7T / dfiKi(X - fi)e r+ x(fi) (39) 

27T J-A r _i 271 J-A r+1 

(r = 1, . . . , JV - 1) where A = A N = and 

(0) (0) _ sin2fc -2A; = (0) K_ 

The integration kernels are obtained from Eq. (|6|) by setting u — 1. The dressed energy of 
the charged excitations is 

£ c (/c) = — 2 (cos k — cosfco) • (41) 

Similarly, one finds reduced integral equations for the spin-components of the density 

n i f^ r ~ 

P,(A) = -^ J ^S r>1 --J_ Ar d^K 2 (X- i A) Pr ^) 

+— / dpK x {X - /i)p r _!(//) + — / dpK^X - fi)p r+ i(ji), (42) 

Z7T J-A r -i Z7T J-A r+ i 

and of the dressed charge matrix (r, s = 1 , . . . , iV — 1) : 

l r A 



1 r is 

£ sr (A) = 8 sr - — dpK 2 (X - pj^srW 

271 J-A s 

+— / d(iKi(\ - /i)^, r _i(/i) + — / dpK x {X - //)£ sr+ i(/i) (43) 

Z7T J-A s _i Z7T J-A.+i 



1 /-A s _i _ X /" A s+1 

dfiK^X - /i)£, r _i(jt) + — / 

V s _i Z7T J-A s + 1 

The other elements of the dressed charge matrix are found to be 

Z cc = l, Z cr = 0, r = 1, . . . , N - 1 



1 /-a 
2^ 



/ Al dA£i(A)£ri(A). (44) 

J— Ai 



The expression for Z rc can be rewritten using the symmetry of the kernel in (O) to obtain a 
simple relation to the densities of electrons with SU(iV) index s: 

1 rAr r n 

Z rc = - dX Pr (X) = l-J2- (45) 



Tl n J~A r s=l 



(since n s = n c /N for h s = this reproduces the corresponding entries in Eq. (|28|) in the u — > oo 
limit). 

For small magnetic fields (corresponding to large but finite values of the A r ) one can employ 
the Wiener-Hopf method to the integral equations ( |39| ) for the dressed energies together with 
condition (|l|) to compute the field dependence of the A r (gj(u) is given in Eq. (|62[) ) : 



7rr /— 7rA, r \ / simrrs/N sinirst/N \ ef' 

sm — exp — J - y s i n7ls /2Ng-(-in/N) ) rfT e ^ 

An analogous computation yields the actual field dependence of the Z rc in Eq. (|^). For /i s <C h c 
we find 

/V — r N N ~ 1 h 

(47) 

where s< (s>) is the smaller (greater) of the integers r, s. 

From (|11|) the bare energies are known to be proportional to the applied magnetic 
fields, hence the A r in fl46|) show the logarithmic dependence on the fields found previously in 
the isotropic Heisenberg spin chain ^8f and the SU(2) Hubbard model |§: 




Ar ~ In ^ (48) 



where h is the typical strength of the fields h r . As is known from the SU(2) Hubbard model this 
strong dependence on small applied fields shows up in the field dependence of the anomalous 
dimensions since they contain terms oc 1 /A s as the leading corrections in their spin-components 
Af (r = 1, . . . , N — 1) (see Eq. fl30D ) and is related to a nonanalytic field dependence of the 
magnetic susceptibility in this system |18] and the SU(iV) Heisenberg model [2J|. This is in 



contrast to the charge-components of the conformal dimensions where fl47D implies a linear 
dependence on the applied fields. 

For sufficiently large fields (h ~ o(h c )) the system saturates in a state with all electrons 
occupying states in the band(s) with the lowest magnetic energy. This final state depends on 
the particular choice of the magnetic fields h s in ([!]). In the following we shall consider two 
possible cases explicitely, generalization to others is straightforward. 

—One natural interpretation of the SU(iV)-index is that of an orbital quantum number, i.e. 
taking the electrons in the r-th band as having spin S + 1 — r where N = 2S + 1. In this picture 
the coupling to the magnetic fields should be through Zeeman terms giving 



for the fields in ([!]) or 



h r = -(S+l-r)h 7 l<r<N = 2S+l (49) 



e 



^ = h - ^y,,, (50) 



for the bare energies (|40"D . In this interpretation it is straightforward to see that the quantity 
h c introduced in Eq. fl4"0"D is simply the large-w limit of the critical magnetic field beyond which 
the ground state of the system is ferromagnetically ordered (i.e. all electrons are in the band 
with spin-S 1 ). For h > h c only excitations with AM r = (r = 1 . . . , N — 1) are gapless. By 
construction, the corresponding correlation functions are those of free spinless electrons. 

—Another possible interpretation of the N bands in the degenerate Hubbard model is that 
of degenerate bands of spin-±| electrons. To be specific let us choose the electrons in the first 



iV + bands as having spin f and the ones in the remaining AL = N — N+ bands as carrying spin 
|. This choice gives 

h h 
K = -- for < r < JV+, h r = - for N + < r < N, (51) 

for the coupling of a physical magnetic field to the system corresponding to 

e (o) = h8 r>N+ - -^^,1. (52) 

1 + A 

for the bare energies (0). Again the system saturates at large fields: for 

fc ^ = ^W5 + ^)-*G)) (53) 

the groundstate of the system is determined by filled bands for the f-electrons with densities 
n s = n c/N + while the N_ bands of j-electrons are empty. This state shows SU(iV + ) spin- 
symmetry, excitations involving creation of electrons in one of the j-bands are massive. The 
gapless excitations and corresponding critical exponents are given by the u — *> oo-limits of 
the expressions in the preceeding section: the dressed charge matrix for the critical degrees of 
freedom is of the form ( pgj) with N replaced by iV+ and z c — 1. 



A Appendix 

In this appendix we briefly list some mathematical results which are helpful in the solution 
of the Wiener-Hopf equations in Sections 4 and 5. The Wiener- Hopf method itself has been 
reviewed in the appendix of 0. 

In the reduction of the matrix Wiener-Hopf problem to scalar ones one has to diagonalize 
a tridiagonal (N — 1) x (TV — 1) Toeplitz-matrix of the type 





(2y 


—x 


•• 


• 0\ 




—x 


2y 


—x 


. 


T = 





—x 


2y ■• 


. 



V o 



(54) 



—x 2y ) 



where x and y are real numbers. In terms of A± = y ± \/y 2 — x 2 one finds for the determinant 
of T 

det T = x + + _ x , (55) 
where we assumed \y\ ^ \x\. Using this result the eigenvalues tj of T are easily found to be 



(56) 



with corresponding eigenvectors 



\Vn-iJ 



2 . fjl% 



(57) 



Thus T is diagonalized by the matrix U G 0(iV — 1) with elements 



^/^(^V (58) 



The inverse of T can be written as follows: 



T' 1 = -r-= S (59) 
det T K J 

where the elements of the matrix S are given by 

A"* — A- 7 

S j i = a j< a N - j> x i> - : ' < , a,j = A ^ _ A ~ ■ (60) 

Here again j < = min (j, I) and j > = max (j, I). 

The main step of the Wiener- Hopf procedure is the factorization of the Fourier-transformed 
kernel into a product of two functions g^ 1 which are analytic in the upper and lower complex 
cu-plane, respectively. In the derivation of ( f46"|) we have used 

S7T 



2e"H (roshcu - cos— j = g+(u)g-{w) (61) 



with 



where T denotes the Gamma-function. Note that the asymptotic behavior of the functions 
defined in (|62|) is simply lim^oo gf{ui) = 1. 



■ exp z^(l + ln2) 

L \2N 2-k) \ 2N 2-k j 



Finally we want note that in the derivation of (|29"D we do not need the explicit form of the 
corresponding kernel K(uj) into factors G ± (co | ). One only needs the value G ± {uj = 0) which is - 



due to the symmetry property G + (lu) = G (—00) - equal to yG + (0)G (0) = JK(0) (see also 

!)• 



References 

[1] F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980); J. Phys. C14, 2589 (1981). 

[2] A. A. Belavin, A. M. Polyakov and A. B. Zamolodchikov, Nucl. Phys. B241, 333 (1984). 

[3] J. L. Cardy, Nucl. Phys. B270 [FS16], 186 (1986). 

[4] H. W. J. Blote, J. L. Cardy and M. P. Nightingale, Phys. Rev. Lett. 56, 742 (1986); I. 
Affleck, ibid. 56, 746 (1986). 

[5] A. G. Izergin, V. E. Korepin, and N. Yu. Reshetikhin, J. Phys. A22,2615 (1989). 

[6] H. Frahm and N.-C. Yu, J. Phys. A23, 2115 (1990). 

[7] F. Woynarovich, J. Phys. A22, 4243 (1989). 



[8] H. Frahm and V. E. Korepin, Phys. Rev. B42, 10553 (1990). 

[9] H. Frahm and V. E. Korepin, Phys. Rev. B43, 5653 (1991). 

[10] N. Kawakami and S.-K. Yang, J. Phys.: Condens. Matt. 3, 5983 (1991). 

[11] M. Ogata, M. U. Luchini, S. Sorella, and F. F. Assaad, Phys. Rev. Lett. 66, 2388 (1991). 

[12] C. N. Yang, Phys. Rev. Lett. 19, 1312 (1967). 

[13] F. D. M. Haldane, Phys. Lett. 80A, 281 (1980); ibid. 81A, 545 (1981) erratum. 

[14] P. Schlottmann, Phys. Rev. B43, 3101 (1991). 

[15] A. Schadschneider, Dissertation, Universitat zu Koln (1991). 

[16] T. C. Choy, Phys. Lett. 80A, 49 (1980). 

[17] T. C. Choy and F. D. M. Haldane, Phys. Lett. 90A, 83 (1982). 

[18] K. Lee and P. Schlottmann, Phys. Rev. Lett. 63, 2299 (1989). 

[19] K. Lee and P. Schlottmann, Physica B163, 398 (1990). 

[20] P. Schlottmann, Phys. Rev. B45, 5784 (1992). 

[21] E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 

[22] A. Kliimper, A. Schadschneider and J. Zittartz, Z. Phys. B 78, 99 (1990). 

[23] A. Schadschneider and J. Zittartz, Z. Phys. B 82, 387 (1991). 

[24] B. Sutherland, Phys. Rev. B12, 3795 (1975). 

[25] C. N. Yang and C. P. Yang, Phys. Rev. 150, 327 (1966). 

[26] H. J. de Vega, J. Phys. A21, L1089 (1988). 

[27] J. Suzuki, J. Phys. A21, L1175 (1988). 

[28] N. M. Bogoliubov, A. G. Izergin and V. E. Korepin, Nucl. Phys. B275 [FS17], 687 (1986). 

[29] P. Schlottmann, Phys. Rev. B45, 5293 (1992). 



Figure Captions 

Figure 1: 



Dependence of the critical value u c the degeneracy iV as obtained from fll3|) . Note that for 
u > u c the system is in an insulating phase whereas for u < u c it shows metallic behavior (for 
n r = 1). 



Figure 2: 



(a) Lines of constant z c for vanishing magnetic fields in the n c -w-plane for N = 2. The drawn 
lines correspond to z c = 1.0904, 1.1679, 1.2237, 1.2622, 1.2892, 1.3088, 1.3234 (corresponding 
to z = 0.4, 0.8, . . . , 2.8). Note that z c — > 1 for u — > and z c — > x/iV for -u — > oo with 
n c arbitrary, n c — > with u arbitrary, and n c — > 1 with w > u c . (b) As in (a) but for 
N = 4, z c = 1.1241, 1.2445, 1.3473, 1.4318, 1.5013, 1.5586, 1.6062. (c) As in (a) but for N = 6, 
z c = 1.1316,1.2630,1.3799,1.4812,1.5689,1.6453,1.7122. (d) As in (a) but for N = 8, z c = 
1.1344, 1.2701, 1.3928, 1.5012, 1.5973, 1.6830, 1.7601. 



Figure 3: 



